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Abstract 

We solve for the S l O(3)-invariant Kahler-Einstein metric on P 2 with cone singularities along 
a smooth conic curve using numerical approach. The numerical results show the sharp range of 
angles ( (n/4, 2n]) for the solvability of equations, and the right limit metric space (P(l, 1,4)). 
" These result exactly match our theoretical conclusion. 

1 Introdution 

Let D be a smooth conic curve in P 2 . In this work, we fix D = {Zf + Z\ + Z\ = 0}. In the recent 
work [4] , we have considered the problem of existence of Kahler-Einstein metrics on P 2 with cone 
'— '< singularities along D of cone angle 2irf3 € (0, 2tt]. The following is the main result in this study [4]: 

£NJ ■ Theorem 1.1 ([4]). There exists a conical Kahler-Einstein metric on (P 2 , (1 — /3)D) if and only 

if/3e (1/4,1]. 

, As pointed out to us by Dr. H-J. Hein, when j3 = |, this gives rise to cone Calabi-Yau cone 

■ metric on the 3-dimensional A2 singularity x\ + x\ + x\ + x\ =0. 

\D : This is a question raised by Gauntlett-Martelli-Sparks-Yau in [3]. In [3], they proved there can 

■ not exist such Calabi-Yau cone metric on 3-dimensional Ak-i singularities xf + x% + x% + x% = if 
k > 4. The idea is to look at the links Lk of such singularities. Any such Calabi-Yau cone metric 
would induce a Sasaki-Einstein structure on L^. By further taking quotient by the U(l) action 

' . I , generated by the natural Reeb vector field, we would get an orbifold Kahler-Einstein metric on 

(P 2 , (l — j)D). In [3], the obstruction for k > 4 comes from the Lichnerowics obstruction. In [4] this 
I was explained as (P 2 , (1 — \)D) being not log-K-stablc if k > 4. For k = 1 and k = 2 case, we have 

the standard examples corresponding to the P 2 with Fubini-Study metric and (P 2 , \D) = P 1 x P 1 
with the product metric. These discussion leaves open the existence problem when k = 3. 

The new insight from [4] is that we can put such kind of orbifold Kahler metrics in the more 
broad family of conical Kahler metrics. In our notation /3 = 1/fc. This allows us to give a uniform 
theory which together with an interpolation argument lead us to Theorem 1.1. 

However, as pointed out in [4], such result is in contradiction to the result by Conti in [1], which 
says there is no cone Calabi-Yau cone metric on A2 singularities. His proof is by classifying all the 
cohomogencity one 5-dimcnsional Sasaki-Einstein manifolds. This leaves us wondering which one 
is right. 

We decide to attack this question by returning to the approach in [3] where the equations of 
orbifold Kahler-Einstein metrics on (P 2 , (1 — l/k)D) were written down. Note that because of 
50(3) symmetry, such equation comes from the work in [2]. Moreover, the transformation and 
change of variables introduced in [3] is very useful for dealing with the problem at hand. In this 
way, we get a 2nd order differential equation with appropriate boundary conditions. 

Since we could not integrate the equation for general (3 we will use numerical simulation to 
solve it. This was suggested in [3]. Our goal is to carry out such numerical approach. As it turns 
out, the result is same as we expected. 
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Theorem 1.2. The equations corresponding to 50(3) -invariant Kahler- Einstein metric up on 
(P 2 , (1 — P)D) has a numerical solution if and only if j3 > 1/4. 

As suggested by Dr. Song Sun, we will further verify the conjecture proposed in [4] which 
predicts the limit metric space as (3 goes to the critical value 1/4. Again, the numerical result fits 
well with our expectation. 

Theorem 1.3. As /3 — > 1/4, the metric space (P 2 , cup) converges to the metric space (P(l, 1, 4), Coke) 
where Coke is the induced orbifold Kahler- Einstein metric coming from the standard Fubini- Study 
metric on P 2 by the natural branch cover: P(l, 1, 1) — > P 2 (l, 1,4). 

The precise meaning of the above statement is detailed in Section 4 and Section 5. This gives 
direct support to our result in Theorem 1.1. We hope to discuss some generalizations of this 
example in [4] later. 

The organization of this note is as follows. The first section gives a detailed review of the 
structure of S'0(3)-orbits for P 2 . The second section discusses the equations we want to solve. 
Again, we carefully review the approach in [3] and work out more details. In the third sections, we 
show our first numerical result Theorem 1.2. In the last section, after describing the SC/(2)-orbits 
of P(l, 1,4) we demonstrate our numerical studies which explains Theorem 1.3. 



2 S0{3) orbits 

Let us first review how to decompose P 2 into 5*0(3, R)-orbits following [3]. First note that P 2 = 
(C 3 — {0})/C* under the equivalence relation (Zi, Z 3 ) ~ (XZi, XZ 2 , XZ 3 ) for some A^OeC*. 
Now fix any / Z := (^i)i=i <= C 3 , ^ determines a point in P 2 with homogeneous coordinate 
[Z] := [Zi]f =1 = [Z\, Z2, Z3]. Now write the polar decomposition 



So if we define 

then [Z ?; ] 3 =1 = [Z t f l=l and 

Now write 



Z( + Z% + Zi = p 2 e 2l °. 

Zi = e~ ie Zi, 

Z\ + Z\ + Z\ = p 2 > 0. (1) 



Z t = m + v^ivi 
then the identity (1) is equivalent to the identity 



M 2 -M 2 =r 2 ; 



u-v = 0. (2) 



We use these two relations to define the set: 

O = {{u,v) ^0\u-v = 0, M 2 - H 2 > 0} C (K 3 ) 2 - {0}. 

Define a equivalence relation on O by (u, v) ~ a(u, v) for some nel*. Denote the quotient set by 
O = O/ ~. Then we have defined a homcomorphism 

$ : P 2 — ► O 

[Zi]i=i ^ [«, v] satisfying u + V^Tw = e -iAr g (z?+z|+z?) ^ z 2 ,Z 3 ). 
The 50(3) acts on P 2 = O by 

9 ■ (u,v) = (gu,gv). 
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The quotient of this action is an interval: 

R : O — > [0,1] 

[u,v] » M 
\u\ 

So the function R classifies 50(3) orbit. Moreover it's easy to verify that equivalently we have the 
relation 

\Z? + Zl + Zl\ l-R* 
\Z^ + \Z^ + \Z^ l + R?- w 

Now at each point (v,w) G O, we can get a orthonormal basis (e u = u/|u|,e„ = f/|f|,e„, := 
e u x et,). We will denote U(l)i, U (1)2 and U(l) 3 to be the rotation around the axes in the direction 
e u , e v and e w respectively. 

Lemma 2.1. The generic orbit is Orb^ = i? = SO(2i)/'Ea (when < Ro — R([u,v]) < I). The two 
special orbits are 

Orb fl=0 = (S , 0(3)/Z 2 )/t/(l)i = MP 2 
Orb fl=1 = (SO(3)/Z 2 )/U(l) 3 = P 1 

Proof. When < R = M < 1, the stabilizer of 5*0(3) action at [v,w] is isomorphic to Z 2 with 
generator being the rotation around e w with angle n, i.e. (e 

When i? = 0, i>=0. The stabilizer is generated by Z2 and U(l)\. The generator of Z2 can be 
chosen to be (e U) e 2l e 3 ) i-> (— e u , — e 2l e 3 ) for any e 2 , e 3 such that {e„, 62,63} is an orthonormal 
basis. U(l)\ is the rotation group around e u . It's easy to verify that 

Orh R=0 = (M 3 - {0})/R* = KP 2 . 

When R — 1, |w| = |u|. The stabilizer is [/(l)-rotation group around e w denoted as 1/(1)3. Note 
Z 2 C t/(l)3. It's easy to see that (for example by (3)) 



Orb« = i = {Z{ + Zi + Zi = 0} = P 1 C 

Fix the generator of so(3) = Lic(SO(3)) to be 
/ \ / 

x x = -1 \ ,x 2 = { I , x 3 



□ 






Then the corresponding invariant vector field on the orbit SO(3)([u 7 v]) at point [u,v] is given by 
the infinitesimal rotation around three axes in the directions of e u , e Vl e w respectively. In other 
words, they are the generator of the action of E/(l)i, U(l) 2l U(l) 3 respectively. 



1. Around e u : 

T = — 

u ~ cW 



1 — l(cos 9e v + sin9e w )\v\) = s/—l\v\ 



8=0 



2. Around e v : T v = ^ I (— sin^e.^ + cos9e u )\u\ + yf^Iv = — \u\e w . 



dfi 1 61=0 

3. Around e w : 

d 



T 

' d6 



(|u|(cos6le u + sin 6^) + y/— 1(— sin#e u + cos9e v )\v\) 

=0 

u\e v - \/^l\v\e u . 
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Tr= To 



( {4 + 



We can define another vector field generating the radial transformation 

d 

i 

9=0 

Note that strictly speaking, the above vectors represent the tangent vector in 

T [u+lv f 2 = Rom(C(u,v),C n+1 /C(u,v)). 

Lemma 2.2. On Orb fl=0 = MP 2 , T u = 0; On Orb R=1 = P 1 , T w = 0. 
Proof. When R = 0, M = 0, so T„ = on Orb fl=0 i?P 2 . When i? = 1, 

T w = |u| n - v^T|«| n = -V=T(« + >/=!«) 

so T w = —yf—l(u + y/—lv) € C • so T w \r=\ = 0, i.e. T w vanishes on the special orbits 

Orb fl= i = P 1 . . □ 

Note that this Lemma also follows from Lemma 2.1 by the fact that U(\)\ is the stabilizer 
group on Orbij = o generated by T u , while ^(1)3 is the stabilizer group on Orb^-i generated by 
T 

3 Equations for 50(3) invariant Kahler-Einstein 

For special metrics g on P 2 , we have the following 

Lemma 3.1. 1. For any Kahler metric g, we have \T u \ g < |I\,| g . The equality holds only on 
the special orbit Orb^ = i = P 1 . 

2. For any SO(3) invariant metric g, \T v \ g = \T w \ g on the special orbit Orb^ = o = KP 2 . 

Proof. 1. Because Kahler metric is compatible with complex structure J = i-, so 

<^ = lfl e ^ = IH e -l g= 14 = jR <i. (4 ) 

\T v \g \\u\e w \g \\u\e w \g \u\ 

2. On the special orbit Orbij = o = KP 2 , v = 0. Let 71 (9) = \u\(—sm8e w + cos6e u ) and 
12(0) = |u|(cos0e u +sin0e„). Then T v = j[(0) and T w — 72(0). Because there exist rotations 
g{9) in SO(3) such that g{9) ■71(6') = 72 (0). So the conclusion follows from invariance of the 
metric under SO(3). 

□ 

Now choose the dual basis of T u , T v , T w } to be one forms given by {dR, o~i, 0-2, 03}. For any 
50(3) invariant Kahler metric on P 2 , {T u ,T v ,T w ,Tfi} is orthogonal. The metric can be written 
in the form 

g = (dtf + a 2 cr 2 + b 2 a\ + c 2 a\ . (5) 

where 

dt = -\T R \ g dR, a = \T u \ g , b=\T v \ g , c=\T w \ g . 

The minus sign in the first identity is to make the special orbit P 1 to sit in the distance location. 
By Lemma 2.2 and Lemma 3.1, we know that 

Corollary 3.1. 1. a<b on P 2 . On Orb R=1 = F 1 , c = 0, a = b. 

2. On Orb fl=0 = MP 2 , a = 0, b = c. 
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Example 3.1. When j3 = 1, then the SO(3) invariant metric is the standard Fubini-Study metric 
on P 2 . We can write it in the form of (5). One way to do this is to recall the following description 
of Study-Fubini metric. Let ^(t) := [Zi(t), Z<i(t), Z^(t)] be a curve in P 2 with the tangent vector is 
f(0) = ((2- 1 (0),^ 2 (0),^ 3 (0)) i y (Z£(0),Z£(0),Z£(0))) G Hom(O p2 (l),C" +1 /O p2 (l)). The length 
ofj'(0) is given by 

wm i M! _ ( l7 , m 2 i(^(o),z(o))i 2 \ 

17 (o)\fs = -j^r - [\z (0)1 — ) /\z(0)\ ■ 



Using this formula, it 's easy to verify that 



\Tr\fs - \ u \i u l\ v [i = tfw, \T u \fs 



\Tv\fS = L,|l"L,|2 = ./-,! p5 ; \ T w\F 



H = R 
VM 2 + M 2 VTTi! 2 
M 2 -M 2 1-i? 2 



M 2 +M 2 v / TW' i- 1 """ |u| 2 + |i;| a 1 + i? 2 ' 
So the normal distance function t is determined by 



alt = & t(l) = => R = tan ( ^ - * 

1 + i? 2 w \4 



5*0 < i < tt/4 and 



i ( — — tj = cos (t + — 



b = sin (t + , c = cos ^ - 2t) = sin(2i). 



Example 3.2. Without further ado, we give the data corresponding to P 1 x P 1 = (P 2 , \D). (See 
[2] and [3]) 

a(t) = cos(^), b(t) = -j=, c(t) = — sin(V3t). 

The range for t is < t < 7r/(2\/3). 

By [2] and [3], the Kahlcr-Einstcin equation is reduced to a system of ODEs: 

" — 2bc 

6 = (6) 

Note that the equation in [3] defers from [2] by a (negative) factor (—abc) which is caused by a 
change of variable. 

The boundary condition is given at t = corresponds to the special orbit Orb_R = i = P 2 where 
by Corollary 3.1 a = \T u \ g = \T v \ g = b and c = \T w \ g = 0. Moreover, the cone angle equal to 2tt/3 
along Orbt = o = P 1 requires c = 2{3. The factor 2 comes from the fact that when < R < 1 the 
stabilizer is Z2. So the boundary is given 

a{t) = a + 0(t) 
b(t) = a + 0{t) 
c(t) = 20t + 0{t 2 ) 

Since the normalized Kahlcr-Einstein metric satisfies 

Ric(oj' ) = ^ + (1 - /?){£}. 
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Because [D] = |ci(P 2 ), so, by taking cohomological classes on both sides, we get 



m = -(l + 2/?). Cl (P 2 ). 

So a and 8 arc related by a 2 = 5 ■ | (1 + 2/3) since both sides are proportional to the volume of P . 
The factor 6 can be carefully tracked out, but it can also be easily determined either by checking 
the standard P 2 with Fubini-Study metric in Example 3.1 or by substituting in to the last equation 
in (6). The result is 

o 1 



a 



6 



(1 



obtained from equation (6). 

When t = t*, we know from Corollary 3.1 that a(t„) = and i>(i*) = c(t*). 



Lemma 3.2. 



o(i*) = -l, b(U) = c{U) = 0. 



(7) 



Proof. From the first equation in (6) and b(t*) = c(i«), we get d(t*) = — 1. Then we use this to 
derive from Equation (6) that 



b(i*) = — c(i*) 
So the 2nd identity follows. 



lim 



-(&>*) - c(U)) = -2b{U). 



□ 



Note that d(t*) = — 1 is compatible with the fact that the metric is smooth along Orb^-o — 
RP 2 . 

Note the solutions of equation (6) is not unique around the point (a(0), 6(0), c(0)) = (a, a, 0). 
There are at least three possibilities: a < b, a = b, a > b. The a = b case corresponds to the 
Gibbons-Pope-Pederson metric as pointed out in [2]. We are in the a < b case. The symmetry of 
a, b is broken by writing down the differential equation for the variable R = a/b. Using (6), we get 



dt 



1. 



So it's natural to do the following change of variables introduced by [3]. 

dr 

Then 



1/c. 



dR 
dr 



— = R Z -I. 



Using a < b ((4)), we get the solution 

R = - = -tanh(r). 

Moreover, we get the range for r: — oo < r < 0. We list the the ranges of R,t,r as follows: 





P 2 


SO(3)/Z 2 


RP' 2 


R 


R = l 


1 > R > 


R = 


t 


t = 


< t < u 


t t j|s 


r 


r = — oo 


—go < r < 


r = 



Define / = ab, then / satisfies the second order differential equation 

(fir ) =2[6/ + coth(2r)]. 



dt 



dr 
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Example 3.3. By easy calculations, one can get that, for P 2 , / = — 4tanh(2r), / r (0) = 1; and 
for P 1 x P 1 , / = -itanh(r), / r (0) = ±. See [3]. 

Let h = f r then this is equivalent to a system: 

fr = h 

(8) 

h r = 12//i + 2coth(2r)ft- ^. 
It's easy to verify that the data (/, R, h) and (a, b, c) determine each other by the relation 

f = ab, R=^, h = f r = -c 2 . (9) 
The boundary condition is given by 

/(-co) = a 2 ,/(0) = 0. 

h(0) = f r (0) = -c(U) 2 = -b(Q 2 . 
Using (7), (6) and t r (t) = c(t), we get 

MO) = /rr(0) = (fttU + ft{tr)t)U) |t=t» = a{U)b{Q 3 = 0. 



4 Numerical Studies: /3 > 1/4 

Now we explain our numerical simulation. We introduce the variable r for convenience and choose 
boundary value (/(0), h(0)) = (0, — - := — &(t*) 2 ) and solve the equation (8) numerically. However, 
this can not be done because there is a zero on the denominator for r = on the second equation 
in (8) (although it's canceled by zero on the numerator). We can however move away from r = 
a little bit by using the boundary condition and Taylor expansion: 

f(r) = f(Q) + fr(0)r + O(r 2 ) = --r + O(r 2 ) 

T 

h(r) = h{0) + h r {Q)r + 0{r 2 ) = -- + 0{r 2 ) 

T 

So numerically, we can choose ro < to be very close to and choose the boundary condition to 
be 

(/(roUM) = (--,--)• 

T T 

For example, in the following numerical simulation, we choose ro = — 10 -5 . Then we can shoot the 
trajectory out for r going from ro backward to — oo. The following numerical solutions correspond 
to P 2 when r = 1 and P 1 x P 1 = (P 2 , ^D) when r = 3 respectively. They can be obtained for 
example by the ND Solve tool in Mathematica. 
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Of course, the above graphs of / = f(r) just recover the graph f(r) = — i tanh(2r) for P 2 and 
f(r) = — | tanh(r) for P 1 x P 1 (up to high precision). 

If we choose different r, then we get different solution f,h — f r . We know that lim r _ i ._ 00 f(r) = 
/(— oo) = a 2 = 1- g • Numerically, we can just evaluate f(r) for r being sufficiently negative to 
calculate a 2 . Actually, after several tests, one can observe that for fixed r, the graph becomes 
flat very quickly which means f(r) becomes stabilized very quickly as r goes toward — oo. We 
can use Mathematica to calculate (very dense) sequences of data for {r, /(r, r)} where we make 
solution / depend the boundary data r. Then we sample the value of /(r, r) at r = —300. (One 
can certainly choose r to be more negative but the visual effect does not change) The following is 
the numerical result. The two pictures are for short range and long range of r respectively. 




We see immediately that a 2 is a decreasing function of r. More importantly, from the picture, 
we sees that one always has 



a 2 = 



1 + 2$ 



> 0.25 



"4 



and all the (3 > j can be achieved. In particular, when (3 = |, where a 2 = j$ = 0.277777..., one can 
find approximate value of t ~ 0.673 from numerical result. In the picture, we have identified three 
special points: (1,0. 5), (3, 1/3) and (0.673, jg) which corresponds to /3 = 1, | and ^ respectively. 
The corresponding graph of / and g = f r for r = 6.73(0 — |) is the following 




Finally, note that we are only interested when /3 < 1, or equivalently when a 2 < 0.5. However 
the picture suggests we can even pass /? < 1 and solve for conic Kahlcr-Einstcin metric with cone 
angle 2tt(3 > 2n along the conic curve. 
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5 Limit as (5 goes to 1/4 

Wc know that SU(2) acts on P 1 naturally. As pointed out in [4], the following embedding is 
cquivariant with respect to the covering homomorphism cf> : SU(2) —> SO(3, K). 



A : P 1 — > P 2 

[U , Ux] ^ pi + Ul,2iU U u i{Ut - 17?)]. 
Here SU{2) acts on P 2 (l, 1,4) by acting on the first two variables: 

g-[U ,U 1 ,V) = [g-(Uo,U 1 ),V]. 

Note that 

A(P 1 ) = {Zf + Zl + Z\ = 0}. 
Fix generators of SU(2, C) to be standard Pauli matrices: 



-l o y ' \ i o y ' ° \0 -i 

Note that the commutator relation [Yi,!^] = 2Y 3 and cyclicly. So by letting Yi = Yj's satisfy 
[Yi, I2] = Y 3 and cyclicly. Let Y/* be the vector fields on P(l, 1, 4) corresponding to the infinitesimal 
actions of Yi. Then we have 

Lemma 5.1. When we restrict to P 1 , A*Y* = T u , A*Y 2 * = -T v , A*Y 3 * = -T w . 
Proof. A(1,0) = (1,0, i) = u + iv with u = (1,0,0) and v = (0,0, 1). So w = u x v = -(0, 1,0). 
A*Y* = (2(U U + UxUi), 2i{U U 1 + U Ui), 2i(U U - U X U X )) 

1. Y? = §(L7i, -Ob), so 

a* y; = (0, - i (t/ 2 - c/ 2 ), 2z(7 c/ 1 ). 

In particular, F*|(i, ) = 1(0,-1) and A*^*^) = (0,-i,0) = ie w . So A^Y,* = T u . 

2. Y 2 * = i(iC/i,if/ ), so 

A* Y 2 * = (2iU U u -(f/ 2 + Ul), 0). 
In particular, Y 2 *| (1)0) = 0), A*Y 2 *| (li0>i) = (0, -1,0) = e w . So A,Y 2 * = -T v . 

3. Y 3 * = ^(iUo,— iUx), so 

A*YT = (z((7 2 - C/ 2 ), 0, -([/ 2 + Ul)). 
In particular, Y 3 * |n )0 ) = f (1, 0) and A^Yg*^^.^ = (i, 0, —1) = — v + iu. So A*Y 3 * = — T w . 

□ 

We can define a function which classifies the S'C/(2)-orbits 

i?:P(l,l,4) — > [0,+oo) 

'|C/ | 2 + |C/i| 2x 1/2 



[[7 ,J7i,F] 1— > 



|V|V2 



Lemma 5.2. TTie generic orbit when < R < 00 is isomorphic to SU(2)/Zi = SO(3)/Z2- The 
special orbit are 

Orb£ =0 = Pt= [0,0,1], 0rb A=oo =P 1 . 
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Proof. If < R < +00, then [U , U lt V] is the same as {y/^l J U Q , y/^Ut, V], j = 1, 2, 3, 4. So the 
stabilizer is isomorphic to Z4. The cases of special orbits are clear. □ 

Now the SU (2)-invariant Kahler metric has the form 

g = dt 2 + a 2 a\ + b 2 cr 2 + (?a\. 

Similar as the example 3.1 in Section 2, we can calculate the induced orbifold Kahler-Einstein 
metric by the branch covering map: 

P(l,l,l) — > P(l,l,4) 

[Z\,Z2,Z$\ M> [Zi,Z2,Z 3 ]. 

Because the metric is SU(2) invariant, to write down the metric we only need to calculate the 
length of the basic vector fields at the the special point (R, 0, 1) in each S'C/(2)-orbit. 

!• T rWo,i) = (1.0,0), \T k 



2 - ^l(iW) = 3(°.--R»°). a = PA* = 27=^- 

3- *?l(iw) = 3(0,^0), 6^ |y 2 *| 9 = i-^. 

4- ^1(^,0,1) = 1(^,0,0), c = |Y 3 *| 5 = ijjy. 
Again, we can transform to the distance function: 

dt = & R(+oo) = => R = tan(7r/2 - t), < t < ir/2. 

1 + R 



By substituting i? into the expression of a, b and c, we get the data for P(l, 1, 4): 

a = b = - sin - tj = - cos(i). 



c = — sin(7r — 2t) = — sin(2t). 

Note that in this case, a/b = 1. This is very different from the case where /3 > 1/4. For the latter, 
a < b except on the special fibre OrbR = i = P 1 where a = b. Moreover, the boundary condition 
now becomes 

a(t)=b(t) = 1/2 + 0(t 2 ) 
c(t) = i< + 0(i 3 ) 

On the other end where f* = 7r/2, a(7r/2) = b(%/2) = c(n/2) = 0. Geometrically, the special fibre 
Orbfl = o — MP 2 shrinks to a point as /? — > 1/4. If we do the same transformation that dr/dt = 1/c, 
the range of r will becomes (—00, +00) instead of (—00, 0) because c(t*) = 0. 

Next we give the numerical results which show that the metric up converges to the orbifold 
Kahler-Einstein metric on P(l, 1,4). 

First we integrate the identity dr/dt = 1/c numerically and plot the relation between the 
boundary value b(t*) 2 = 1/r and i max = i*. We see that the maximal value for t is an increasing 
function of r. As r — > +00, or cquivalently as ft — > 1/4, t max = t„ converges to 7r/2. 
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Note that the coordinate t is the distance function from the special orbit P 1 . So t is a geometrically 
meaningful coordinate in contrast with r which is only an auxiliary coordinate. So we can get a 
good convergence only when we look the data as functions t. 

Now we can plot the graph of the data set (/ = ab,R = a/b, —f r = c 2 ) as the function of t 
instead of r. (See (9)). The following picture shows the data for four r's: r = 10' for i = 1, 2, 3, 4. 
The corresponding colors and markers are "Blue Round" , "Green Square" , "Orange Triangle" , 
"Pink Diamond" for i — 1,2,3,4 respectively. The "Red Upside-down Triangle" represent the 
data for P(l, 1, 4) where 



f(t)=a{t)b{t) = \ cos 2 (t), 



c 2 (i) = -sin 2 (2i). 



4 v 




0.12 
0.10 - 
0.08 - 
0.06 - 
0.04 
0.02 



One can see that the data for t large fits with the data for P(l, 1,4) very well. Again, we know that 
t going to +oo is equivalent to /? going to 1/4. So the numerical result implies the expected result: 
as (3 — > 1/4, the metric up converges to the orbifold Kahlcr-Einstein metric Cjke on P(l, 1,4). 
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